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ABSTRACT 

Microlensing of extended stellar sources in the LMC and Galaxy by low mass lenses can 
produce variable polarisation. The characteristics of the polarisation and flux profiles 
can provide considerable information about the lens geometry, and help determine the 
importance of such low mass objects as a dark matter component. 
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Recently the search for gravitational microlensing events has 
resulted in a number of successful identifications of lensing 
events of stars in the LMC and the Galactic Bulge. In the 
^ case of the first event reported by the MACHOs group, it 
I would seem that a low mass and possibly compact object in 
the Galactic Halo was responsible for the microlensing (Al- 
lf} cock et al. 1994). The determination of the mass and spatial 
Q\ density distribution of such objects would have important 
**^ implications for assessing their relative importance as a dark 
m atter component (Paczynski 1986). The principal effect of 
' | lensing is amplification of the flux, which essentially depend- 
O s on the distance in units of the Einstein radius of the star 
(source) from the lensing object projected onto the source 
^9 plane. Various methods have been suggested for determin- 
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ing the Einstein radius of the lens, the most promising being 
parallax (Gould 1994a), measurement of spectral shifts for 
rotating stars (Maoz and Gould 1994), and analysis of the 
light curve (Nemiroff and Wickramasinghe 1994). The lat- 
ter two effects depend on the star being considered as an 
extended source. (Even when the radius of the star is con- 
siderably less than the Einstein radius these effects can be 
important). Since the Einstein radius for a given source dis- 
tance depends on the distance and mass of the lens it is im- 
possible to determine both the lens distance and mass from 
the light curve. Of course probabilistic, that is maximum 
likelihood, arguments can be used, although any inferences 
made in this way will be subject to large uncertainties. 

If the source star diameter is larger than a few Einstein 
radii projected onto the source plane, there will be only a s- 
mall variation in the flux amplitude during the lensing event. 
This will make detection of low mass lenses difficult. The 
fact that one has to consider extended sources could poten- 
tially provide more information about the lensing object, 
and could under certain assumptions such as the star's ra- 
dius and the transit velocity of the lens, yield both the mass 



and the distance of the lens (Nemiroff (1994), Witt and Mao 
(1994), Gould (1994b)). This could be achieved not only by 
the analysis of the light curve, but more interestingly, by 
studying the variable polarisation induced by the gravita- 
tional lensing (Simmons et al. 1994), and in this paper we 
address this question. 

The probability that any given star is lensed is ex- 
tremely small. The surveys being carried out depend on 
monitoring millions of stars, and selecting candidate event- 
s on the grounds of a number of criteria, the obvious ones 
being that the light curve is not periodic, and has the expec- 
ted profile of a lensing event. Potentially the statistics of the 
lensing events could yield a lot more information about the 
spatial and mass distribution of lensing objects, although the 
interpretation of these statistics will depend on the validity 
of the assumptions made about the pointlike nature of the 
source star. 

There are three main effects of the star being an exten- 
ded source that we shall discuss in this paper. The first has 
been considered recently in relation to MACHOs by Gould 
(1994a), and by Nemiroff (1994), Witt and Mao (1994). 

(i) There will be a modification of the amplification func- 
tion and thus of the light curve of a lensing event. Our c- 
alculation show that for a source star of radius R > 2.5r]o, 
where rjo is the Einstein radius projected onto the source 
plane (see fig 3) the amplification is always less than 1.34 
(the amplification for a point source at one Einstein radius), 
no matter how close the lens passes by the centre of the star. 
This could be crucial for detecting low mass lenses. Indeed 
it maybe impossible to detect these. 

On the other hand for R < d, where do is the impact para- 
meter on the source plane, the maximum amplification dur- 
ing an event is greater than for a point source. This might 
explain the light curve anomaly for the first event repor- 
ted by the MACHOs group, which shows a high maximum 
amplification (see also Witt and Mao (1994)). 
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(ii) The gravitational lensing of extended sources could 
induce a variable polarisation during the lensing event. This 
was first discussed by Schneider and Wagoner (1987) in re- 
lation to supernovae, and recently by (Simmons et al. 1994) 
in relation to microlensing surveys. The mechanism is the 
following. Light emerging from the limb of the star will nor- 
mally be polarised, as was first discussed by Chandrasekhar 
(1960). During the microlensing of the star there will be an 
amplification of the light from the limb of the star nearest 
to the projected position of the lens, and hence the integ- 
rated light from the source should display a net polarisa- 
tion. Furthermore the position angle of the polarised light 
would simply depend on the angle made by the line join- 
ing the source and lens centres projected onto the source 
plane. Chandrasekhar's calculations for an electron atmo- 
sphere gave a limb polarisation of about 10%. More accur- 
ate calculations indicate a lower degree of polarisation in the 
visible band for grey atmosphere, although high polarisation 
could arise in stellar envelopes (Collins and Buergerl974). 
The maximum degree of polarisation attained during a lens- 
ing event would depend on the star's radius and the distance 
do, or impact parameter, but would typically be less than 
than 10%. Successful measurement of this variable polarisa- 
tion yields still more information about the lensing geometry 
and provides a further means of determining the mass and 
distance of the lensing object. 

One obvious problem with this method is the degree of po- 
larisation and brightness of the star required for the effect 
to be observable. It is however an interesting feature of this 
approach that the rise in polarisation will usually take place 
later than rise in total flux. Typically, since the polarisa- 
tion variability is a differential effect, somewhat akin to tidal 
forces, the period over which the variable polarisation takes 
place will be half that of the flux variability. Thus it would 
be possible to implement an early warning system: use the 
detection of flux variability to alert a search for polarisation 
variability on a larger telescope. 

(iii) Finally, although for point sources and for extended 
sources with no limb darkening the light curves should be 
independent of the wavelength of light, or achromatic, it 
is possible that chromatism could come about through the 
wavelength dependency of limb darkening. For most situ- 
ations we would not expect this to be an important effect, 
however. 

With the detection of more events, the interpretation 
of their statistics will provide more information about the 
distribution of lensing objects. With so few events confirmed 
a full analysis is probably premature, and would require a 
full understanding of the selection effects at play. In this pa- 
per we shall only make a few remarks about the fraction of 
events that would be expected to show polarisation variab- 
ility, and leave a fuller analysis to a future publication. 

Outline of paper is as follows. In section 2 we devel- 
op the theoretical formulation of the problem and intro- 
duce our notation. Our numerical results for polarisation 
and flux variability are presented in section 3. In section 3.4 
we discuss the possible chromatism introduced by extended 
sources and give numerical results for the grey atmosphere, 
and in section 3.5 we make a few remarks about the statistics 
of the events. 
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Figure 1. Two dimensional representation of lensing geometry. 

2 THEORETICAL FORMULATION 

In this paper we shall only consider Schwarzschild lenses, 
where the lensing object has a spherically symmetric grav- 
itational field. 

2.1 The lens equation 

Two images of any point source on the source plane are 
formed by the lens, except in the case where the source point, 
lens and observer are all aligned, in which case an Einstein 
ring is formed. A light ray, see figure 1, is deflected through 
an angle 4GM/rc 2 or 2Rs/r where M is the mass of the 
lens, Rs is the Schwarzschild radius of the lensing object, 
and r = aj,(9. The apparent directions of these images are 
given by 8\ and 82. It is easily shown that these satisfy the 
simple lens equation 

8 2 - 8(6 - do 2 = (1) 

where fi is the actual direction of the source point (i.e. in 
the absence of lensing) and 80 the so called Einstein radius, 
given by 

8o 2 =(- -) 2R S (2) 

\a L a s J 

where 

Rs = mRsQ (3) 

where m = M/M@ and Rsq is the Schwarzschild radius of 
the Sun. 

2.2 Amplification 

We shall introduce as a natural unit of distance the Einstein 
radius projected onto the source plane, rjo, given by 

Vo =C 1/2 (1-0 1/2 ™ 1/2 D (4) 

where D = tJ2Rsqo,s and £ = aj, / as- 

This quantity rjo is fundamental in this analysis. If the 
distance, d, of the source from the lens projected onto the 
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Figure 2. Transit times in days as a function of the fractional 
distance to the LMC assuming that v_\_ = lOO&ms - 1 for a variety 
of lens masses 




lens 



Figure 3. Three dimensional representation of the lensing geo- 
metry illustrating amplification 



source plane is much greater than rjo , there will be negligible 
amplification. If the actual radius of the source star is much 
less than rjo then it can effectively be taken as a point source. 
The time duration of a lensing event can be approximately 
taken to be 



9 *?° o i 
2 - — = 2 v±£ 

t;v j_ 



-3/2 



(1 



(5) 



where v± is the transverse velocity of the lens (see fig. 2). 

If this time is less than two days or so, there is little 
possibility of observing such an event. This puts constraints 
on the source distance and the lens mass for a given lens 
velocity and source distance (see figure 7). 

Amplification of the source essentially arises from the 
fact that the apparent solid angle, dQ,' subtended at the 
observer by an element of area of the source is greater than 
the solid angle, cifi, subtended by the same element in the 
absence of the lens (see figure 3). The specific intensity, /, 
along the (deflected) ray, however, remains constant. Hence, 
with reference to fig. 3, the flux at the observer due to this 
element is amplified and may be expressed as 

dF = I((,<f>)dQ' 

a s z 



(6) 
(7) 



where we have of course assumed that the angles 9, do, 
and fi are very much less than 1, and A is the amplification. 
Because of the azimuthal symmetry of the lens, this amp- 
lification factor can simply be written as a function of the 
distance of the source element from the projected position 
of the lens in units of jjo, C> an d is given by -^6. Since for 
our case of microlensing 

the time delay between the two images is negligible, 
the total amplification of the source element is given by the 
sum of the amplification for the two images. This gives the 
amplification factor as 




Figure 4. Light emerging at angle x to the normal to the stellar 
surface will be limb darkened and polarised. 



A{z) = 
where 



1 ( 1 
- [z + - 

2 V z 



1 + 



(8) 



(9) 



2.3 Extended source 

We shall assume the star to be spherical (see figure 4). The 
specific intensity of the radiation that emerges from the s- 
tellar photosphere will depend on the cosine of the emergent 
angle, fi = cos^. We shall assume this to take the simple 
form 



I(fi) = i + ii ft 



(10) 



where io and i\ can be related to the luminosity of the star 
and the limb darkening (see Chandrasekhar (I960)). The 
flux at the stellar surface is simply 
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position 



Figure 5. The total polarised flux is obtained by rotating the 
contribution from each source element through angle tt — ip and 
integrating over the source 



2ir / I(ji)ndn 
'o 



2 . \ L 



(ii) 



where L is the luminosity of the star and R its radius. The 
ratio ii/io can be found by modelling stellar atmosphere 
using the theory of radiative transfer, and determines the 
degree of limb darkening. For a grey atmosphere, this ratio 
takes the value |. 

The state of polarisation of this emergent radiation is 
best expressed in terms of the Stokes parameters U and Q 
(Chandrasekhar 1960). Light emerging at angle \ 1° t ne 
normal to the stellar surface will be limb darkened and po- 
larised. 

In the coordinate frame whose z axis lies in the plane 
defined by the normal to the photosphere and the emergent 
direction (see figure 5) the Stokes parameters of the emer- 
gent radiation are given by 



U 

Q 



«o(l - fl) 




(12) 



This simply states (assuming Mo to be positive) that the po- 
larisation is in a direction perpendicular to the plane defined 



by the normal and the emergent direction. The degree of po- 
larisation is given by 



P 



(U 2 +Q 2 ) 
I 



2x1/2 



(13) 



From a somewhat simplified c- 

alculation assuming Thompson scattering, Chandrasekhar 
obtained a value of 11.7% for the maximum polarisation, 
that is for light emerging at the limb. 

To obtain the total flux and polarised flux at the ob- 
server we must integrate the contributions of each surface 
element over the source's disc. For the total flux this simply 
gives 



Fi = -^ / (io +•!/*) A(CKd(dcf> 

a * J ./disc 
It is convenient to write equation (14) as 
1 



F 



I = —^(iooio(R, d) + i 1 a 1 (R,d)) 



where 
a (R, d) 

ai(R, d) 



disc 



A(C)CdCd<i> 



M(C)C<W 



(14) 

(15) 

(16) 
(17) 



disc 



In order to evaluate the polarised flux at the observer, we 
need to evaluate Fu and Fq. When summing the contribu- 
tions of the surface elements of the source, we must refer the 
stokes parameters to the same reference frame. We shall take 
this to have its z axis aligned with the line in the source plane 
connecting the source's centre to the projected position of 
the lens (see figure 5). Using the transformation properties 
of the Stokes parameters we obtain in this new frame 



dUnew = d^old cos 2 ^ ~~ ^Qold s i n 2 ^ 

= rft^old cos 2 ^ 

dQnew = d^old s i n 2 ^ + ^Qold cos 2 ^ 

= rft^old s i n 2 ^ 



(18) 



(19) 



Integrating over the disc to obtain the polarised flux at the 
observer we obtain 



Fu = 

and 

F Q = - 
i 

where 



u (l - fi)A(() cos2ip(d(d<f> 



disc 



;U {fi {R, d)-^{R, d)) 



'-'f 



disc 



1 1/2 



m (1 - fi) A(C) sin 2ipC, d(d<f> 



and a 2 = d 2 + ^ 2 + 2£d cos <j> 



(20) 



(21) 



(22) 



and we have introduced the functions 

(A(() cos 2tl>dC,d4> 



MR, d) 
Pi(R, d) 



disc 



CttA(C) cos 2ipd(d<f> 



(23) 



disc 
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Normalised Flux for a range of source/lens configurations 
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Figure 6. Contours in R - do plane of constant flux amplification 
for i\jio = 3/2 corresponding to the grey atmosphere. Amplific- 
ation is expressed in negative magnitudes. Radius of star, R, and 
impact parameter, do , are both expressed in units of rjo . 



The angle ip can similarly be expressed as a function of 
£ and <f>. (The expression for Fq is necessarily zero). When 
the lens lies outside the disc, the integration is considerably 
simplified by taking the origin to be the centre of the disc 
of the source. We shall numerically evaluate these integrals 
and make no attempt to write down analytic forms for them, 
which would necessarily be very complicated. 

To obtain the time profiles of the flux and degree of 
polarisation, we only need note that the distance, d, is given 
by 



d ■■ 



do + [vi_—t 



-11/2 



(24) 



Here do is the impact parameter, or distance of closest ap- 
proach, and v± the transverse velocity of the lens. 



3 NUMERICAL RESULTS 

3.1 Flux as function of radius of star and distance 
of lens 

The results of the flux integration are given in figure 6 for 
the cases of uniform intensity (i\ jio = 0) and limb dark- 
ening corresponding to the grey atmosphere (ii/io = 3/2). 
The flux has been normalised to the flux in the absence of 
lensing. Flux contours are shown on the R and d plane, both 
expressed in units of rjo • There are a number of interesting 
features. For R = 0, the amplification is evidently that a 
point source. For fixed d the flux increases with increases 
radius, until it reaches a maximum at around R = d. (This 
point was observed by Gould (1994b) and Witt and Mao 
(1994)). This means that the probability of detection is in 
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Figure 7. Observable lensing events are constrained by transit 
times and flux amplification. Graphs show the mass of lens that 
produces observable events as a function of fractional distance 
to the source (assumed to be the LMC). Only events above the 
curves are observable. v±_ is taken to be 100 km s~ 1 . 



fact greater than would have been supposed if we had taken 
the sources to be point like. However, as R increases beyond 
d, the flux falls off. For values of R > 3r]o, the flux amp- 
lification is always less than 10% no matter the value of d. 
Thus such lensing events in which the source radius in great- 
er than 3 would not be detectable according to the criteria 
adopted by the current surveys. This would discriminate a- 
gainst low mass and distant lenses (low rjo) being detected 
(see figure 7). 

The time profile of the flux amplification may be ob- 
tained by considering the same figure, but now allowing d 
to decrease to some minimum value, do, and then increase 
again. These time profiles will be discussed below. If d at- 
tains a value less than R, the lens passes over the face of the 
star, and we shall call such an event a transit. If this does 
not occur, i.e. do > R, we call the event a bypass. 



3.2 Polarisation contours 

Figure 8 shows the contours for the degree of polarisation 
in the R and d plane. Polarisation is always perpendicu- 
lar to the line of centres, ie Fu is positive. Once again 
the maximum polarisation for fixed d is achieved at around 
R = d, which corresponds to the maximum amplification of 
the limb. Because of dilution by the rest of the source, the 
value will not reach the Chandrasekhar value at the limb 
of 11.7%. Maximum polarisation will evidently be achieved 
when R is small but of the same size as d. This requires a 
large value of rjo- Such events would not take place very fre- 
quently. However, values of polarisation of more than 0.1% 
would be relatively frequent, and for bright stars (and big 
telescopes) could be observed. 

To obtain the time profiles of polarisation, keep R fixed 
and allow d to decrease to its minimum value do and then 
increase again. Bypasses will show only a single maximum, 
whereas transits will display a double maximum, roughly 
corresponding to when the lens passes over the limb. All 
profiles will of course be symmetric. 
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3.3 Time profiles of flux and polarisation 

As we have already mentioned, lensing events of extended 
sources can be classed into two types, (i) transits and (ii) by- 
passes. In the former, the lens passes over the source whereas 
for bypasses it does not. Case (i) takes place when R > do 
and case (ii) when R < do. The event can be recorded either 
as (a) variation in flux but no polarisation (b) variation in 
polarisation but no variation in flux (c) variation in both 
flux and polarisation. The profiles depend on 5 parameters 
for a given source distance as- These are the radius of the 
star R and the impact parameter do both expressed in unit- 
s of rjo, and v±as/a,L where v± is measured in units of rjo 
per second. (We could choose not to model the stellar atmo- 
sphere, but instead to determine the parameters io, ii and 
Mo from the data). 

If we assume that the radius of the star can be de- 
termined from its spectral type, then rjo can be obtained in 
principle from the time variation of the flux. If we further 
assume a value for the actual transverse velocity of the lens, 
although there is no reason to do so, then the distance to 
the source and the mass of the lens can be determined. 

The fitting of polarisation data will in addition to these 
parameters yield the direction of the transverse velocity. 
Parameter fitting will be particularly easy in the case of 
transits, since do ~ i?cos$ and tv±as/a,L ~ i?sin$ where 
$ is the position angle at limb crossing. 

Examples of time profiles of the flux variation and po- 
larisation variation are given in figures 9 and 10 for various 
combinations of parameters, and correspond to transits and 
bypasses. We leave the scale for time variation arbitrary. For 
comparison, we show the profiles that would be obtained 
for a point source set-up with the same values of d, M, 
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Figure 9. Flux and polarisation profiles for a lens of O.IMq and 
source radii of 10 and 5 Rq for fractional lens distance, £ = 
a L/ a S> °f 1/11 and 10/11. Sources are taken to be in LMC, and 
the time unit is arbitrary. 



f = ah/as and v±. (The flux profiles are similar to those 
obtained by Witt and Mao (1994)). Purely for illustrative 
purposes, we have taken the values of these parameters for 
the MACHO event quoted in the literature (Alcock et al. 
1994). 

For this event, curve fitting of flux data by Alcock et al. 
(1994) gave a maximum amplification of A = 6.86 ± O.f f , 
implying (from eqn. 8) a value of do = 0.15 (in units of rjo), 
and a duration time defined through eqn. (5) of 33.9 ± 0.26 
days . These authors also inferred a most probable value 
of M = O.f 2 Mq for the mass of lens. Let us take the lens 
mass of Alcock et al. and the value of the impact parameter 
for this event to illustrate the utility of polarisation and flux 
measurements for determining the lens geometry. In fact the 
time profile of the flux variation would be identical for a lens 
of the same mass but at a distance of f — aj, /as, that is for a 
lens much closer to the LMC in this case if the source was a 
point source. For extended sources this will not be the case. 
We take cll/cls to be 1/11 and 10/11 and the radius of the 
source star to be f and 5 Rq. v± is taken to be the same for 
all events. The corresponding profiles are shown in figures 9 

The flux profiles for a nearby lens (£ = cll/cls = 1/11) 
and distant lens (£ = 10/11 are indistinguishable for both 
source radii and effectively the same as for a point source 
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Figure 10. As in figure 9 but for lens of O.OOIMq. Notice that 
for most cases £ = 1/11 and £ = 10/11 are now distinguishable. 



since the Einstein radius is an order of magnitude larger 
than the stellar radii. 

The curve displaying the greatest amplification corres- 
ponds to a nearby bypass. The amplification is measurably 
larger than for a point source (Witt and Mao 1994). The po- 
larisation is negligible for pointlike sources, but shows single 
maxima characteristic of bypasses for the two cases where 
the lens is near the source and rises after a measureable 
change in flux. 

If we take the lens mass to be only O.OOIMq, the s- 
ituation is quite different (see figure 10). Here the whole 
gamut of behaviour is displayed. Transits are now seen as 
flat topped flux profiles, and double maxima polarisation 
profiles. From the polarisation profiles for transits it is im- 
mediately apparent that one source is twice the radius of 
the other. For most of these cases the polarisation should be 
measureable. 

Characteristically the rise in polarisation takes place 
later than the rise in the flux. This could be important in 
setting up an early warning system. Identification of an in- 
crease in flux could be used to signal a search for variable 
polarisation. 




Figure 11. Flux amplification profiles at different wavelengths 
for grey atmosphere for a = hv jkT = 1 and 5. 



we saw in the case of the flux and polarisation. Consider 
the case of the monochromatic flux, F v , and suppose for ex- 
ample that the specific intensity, is a function of the 
cosine, fi, of the emergent angle 

The flux amplification at any given frequency will then 
be given by 



A v (d.R) 



F Vl (d,R) _ f J disc WM0CdCd^ 



F Vl (oo,R) 



J7 d 



(25) 



If the specific intensity had the separable form I v = 
f(v)G(ii) then evidently the amplification given by equa- 
tion (25) would be frequency independent, and no chromat- 
ism would be observed. If on the other hand this is not the 
case, a degree of chromatism can be expected. How large 
this effect is will of course depend on the nature of the stel- 
lar photosphere. In the idealised case of a grey atmosphere, 
where the opacity is taken to be frequency independent an- 
d where LTE is assumed, the solution of the equation of 
radiative transfer 

H-?-I v = I v - B v (T(t)) (26) 
or 

yields 
where 



e- T ^B v (T(r)) 



dr_ 



(27) 



B v (T(r)) is the Planck function and T(r) = 
(1/2) 1/4 T e (l + |) 1/4 . T e is the effective temperature of the 
star. 

The chromatism is not strong, and is only evident for 
small values of d (see figure 11). (However, both the as- 
sumptions that the opacity is frequency independent and 
that LTE obtains are dubious for realistic stars). Thus it 
is unlikely that this source of chromatism would affect the 
selection of lensing events. It might possibly be used as an 
additional technique for determining the lensing parameters. 



3.4 Chromatism 

Any phenomenon in the stellar atmosphere that is limb de- 
pendent will be affected by gravitational microlensing, as 



3.5 Statistics of events 

A simple expression can be obtained for the rate at which 
events occur (e.g. Paczynski (1986)). The number of such 
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surface density is homogeneous, we may write the above in- 
tegral 28 as 



Figure 12. Any stellar source will provide an event during time 
interval (0,T) if it lies within the strip < y < wj_X/£ and 
satisfies selection criterion. 



events that are transits or bypasses, polarised/unpolarised 
for a fixed radius of the source is proportional to the length 
of the corresponding interval in the values of do, or impact 
parameter, at which the different types of events can take 
place. If we assume that all events where the amplification 
of the flux is greater than 1.34 are detected, and similarly 
all events showing a maximum polarisation of more than 
0.1% are recorded as polarisation events, then the fraction 
of recorded events of different types can be obtained from 
figures 6 and 8. If for instance, the source radius is 100 R@, 
and the MACHOs are assumed all to have an Einstein radius 
of 500 Rq, then ~ 40% are polarised. For R > 3r]o, there 
will only be polarisation events, whereas for R << rjo, only 
a small fraction of events will be polarised. The distribution 
of events of different types will thus yield information about 
the spatial and mass distribution of the lensing objects, and 
the size distribution of the stellar sources. 

In fact the predicted rate of recorded flux events for 
extended sources will differ considerably from the predicted 
rate for point sources. No flux event will be recorded if the 
radius of the star is greater than about three times the Ein- 
stein radius, rjo- This effectively acts as a selection on the 
mass of the gravitational lenses that will be observed. 

Consider the case of a gravitational lens at distance 
o-L = £,a.s of mass m. We shall also take the radius of the 
source stars to be fixed. In a time interval, T, assumed long 
compared with the typical lens transit time, the number of 
recorded events will be simply given by 



dy 



dz l£(y , z) S (r]o , R, z,v±) 



(28) 



where £(«/, z) is the surface number density of source stars 
(and \z\ = r/odo). 

S(r]o, R, do, v j_) is the selection function. If we assume 
for simplicity that selection is only on magnification amp- 
lification, then for a point source S = &(r]o — do). For an 
extended source this is no longer the case. Flux amplifica- 
tion will surpass 1.34 roughly when R < arjo and r]o — do < 0. 
(a PS 3) (see figure (12). Thus in this case the selection func- 
tion may be written S PS &(r]o — do)&(r]o<x — R), where 
is the Heaviside function. In the case where we assume the 



T 



v j_ Sso / dz S(r]o, R, z, v j_) 



For a point source this gives an event rate of 
Sso / dz Q(r]o - \z\) = jv± E so Vo 

^ J — oo ^ 



(29) 



(30) 



and for an extended source, with the approximate selection 
function given above, yields an event rate of 

2— E so Vo @(«»?o - R) 
C 

If we now assume that the number density of lenses, raj,, 
depends only on the radial distance to the sources (LMC 
say) then we may write the total rate of events as 

2(v L )a s ( alnsoAQ, j dff n L {£) i]o(0 0(a»?o(O - R) 

= 2{vj_}a s Ns I dff J?o (O™i(O («'?o(O-#) (31) 
Jo 

where (v±) is the mean transverse velocity. Now rjo is given 
by equation (4) as a function of m and If we take the 
number density of lenses to be constant, then this expression 
reduces to 



2(fj_) a s Nstil I ' d{{»7o(f) 



(32) 

where p = R/(amD) is a dimensionless radius of the star 
and D = ^/2Rsq (is as before. This can be integrated ana- 
lytically to yield 

2(v±.)a s Nsn L DF(p) (33) 
where 

A normalised plot of this function is given in figure (13), 
and shows the fraction of lensing events that would actually 
be detected if the star was an extended source compared 
with the number predicted assuming that the source was 
point-like. For source stars of 10 solar radii and a = 3 only 
about 5% of events corresponding to lens masses of 10 -4 
solar masses will be picked up, whereas for lens masses of 
0.1 solar mass about 99% will be detected. 

This analysis could be trivially extended to cover the 
case where there was a distribution of lens masses, and a 
distribution of stellar radii of source stars. 

(Obviously, the rate at which events are detected will 
depend on a number of selection effects, which we shall not 
attempt to model). 



4 CONCLUSIONS 

Microlensing of stellar sources should produce variable po- 
larisation. Depending on the geometry of the lensing sys- 
tem, this could be as high as 2%. A significant fraction of 
MACHO events in which stars from the LMC are lensed 
and manifesting flux variation should also display variable 
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Figure 13. If the sources are extended rather than pointlike then 
fewer lensing events will take place. The fraction of observable 
events for an extended source compared with the predicted num- 
ber of observable events for pointlike source as a function of the 
dimensionless source radius, p = R/ (amD) 



polarisation above 0.1% polarisation, although the precise 
number will depend on the size distribution of stars, and the 
mass and spatial distribution of lenses. Evidently, for nearby 
lenses the fraction will be small, and polarisation events are 
unlikely. For lensing objects lying closer to the LMC, polar- 
isation events and transits will be more frequent. 

In principle the measurement of the variable polarisa- 
tion yields the Einstein radius of the lens and its velocity 
direction on the sky, and thus provides more information 
than simple measurement of the flux variation assuming an 
extended source. Transits would be easily distinguished by 
their double peak polarisation profile. 

The low levels of polarisation produced by lensing mean 
that the stars would have to be sufficiently bright for a giv- 
en size of telescope and given integration time. It should 
be feasible to measure 0.1% polarisation of a star brighter 
than 18 magnitude with one hour's integration time on a 
1 metre telescope. The intrinsically brighter stars are more 
likely to display higher levels of limb polarisation and are 
also more amenable to observation. 

During lensing the rise in polarisation generally takes 
place significantly later than the rise in flux. This provides 
the opportunity of implementing an early warning system 
for detailed study of the polarisation variability. This also 
in part offsets the difficulty in observing low degrees of po- 
larisation for faint stars, since a larger telescope could be 
brought in for the observation of polarisation. Since we are 
here dealing with variable polarisation, the effect of inter- 
stellar polarisation could also be accounted for. 

For high mass lenses variable polarisation would be 
more likely to be observed for more extensive stars. There is 
also the interesting possibility of probing with gravitational 
lenses the stellar atmosphere and envelopes, which can be 
quite extensive, of bright young stars. 

We have presented the simplest model for limb polar- 
isation in stars. More detailed modelling of the stellar pho- 
tosphere for stars of given spectral type should allow more 
accurate determination of the lensing parameters, although 



in principle the parameters determining the limb depend- 
ency of the polarisation could be fitted as well. Thus po- 
larisation measurement should be a valuable supplementary 
tool in the analysis of microlensing by compact objects in 
the galaxy. 
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